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Abstract 



. We examine the theoretical foundations of standard methods for computing 

. density perturbations in inflationary models. We find that: (1) the time-delay 

o . 

. formalism (introduced by Guth and Pi, 1982) is only valid when inflation is 

well-described by the de Sitter solution and the equation-of-state is nearly 
Q . unchanging; and, (2) the horizon-crossing/Bessel approximation extends to 

In: 

c/3 . non-exponential inflation, but only if the equation-of-state is changing slowly. 

Integration of the gauge-invariant perturbation equations mode-by-mode is 
the only method reliable for general models. For models with rapidly varying 
equation-of-state, the correction leads to significantly different predictions for 
the microwave background anisotropy. An important corollary is that meth- 
ods proposed for "reconstruction" of the inflaton potential from anisotropy 
data are unreliable for general models. 
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One of the most important predictions of inflationary cosmologyoLl is that quantum 
fluctuations of the inflaton field grow into cosmological energy-density perturbations.iHi. 
In this paper, we analyze and compare the standard methods for computing perturbation 
spectra in inflation. This consideration is motivated by the need to have theoretical predic- 
tions which match the precision anticipated in forthcoming measurements. We find that the 
simplest, most commonly used methods are approximations with narrow ranges of validity. 
The only reliable method for general potentials is the gauge invariant methodi'i0 in which 
the equation-of-motion for the perturbation must be integrated for each mode. 

The paper has several disparate components which we have organized by section for 
the convenience of the reader: (1) a review of the "exact" gauge invariant methodJffi 
with attention to some subtleties which have caused confusion in past literature; (2) an 
analysis of the time-delay formalism demonstrating that it is valid only when inflation is 
nearly de Sitter (exponential inflation); (3) an analysis of the horizon-crossing/Bessel func- 
tion approach,ii'3ll3 which show that it extends to non-exponential inflation but only if 
the equation-of-state is changing slowly; (4) example (Figure 1) of how the approximate 
methods can lead to large errors in computation of cosmic microwave background (CMB) 
anisotropy in cases of rapidly varying equation-of-state; (5) a summary in terms of conditions 
on the inflaton potential for applying the time-delay and horizon-crossing approximations 
and discussion of implications for "reconstructing" the inflaton potential and constraining 
cosmological parameters from (CMB) anisotropy and large-scale structure data. 
Exact Methodii'Uli In this paper, we consider the case of a single inflaton field; multi- 
field infiation introduces other subtleties that we will discuss in a future paper.El The most 
general form for the metric with linear scalar perturbations is 

rfs2 ^ a2(r){(l + 20)dr2 - 2B\idx'dT - [(1 - 2^)5^^ + 2E\ij\dx'dx^} (1) 

where 0, 5, ip and E are arbitrary functions of space and time. A gauge-invariant combina- 
tion of these variables is the gravitational potential. 



^ = ct> + -^[{B-E')at (2) 

where prime means derivative with respect to conformal time r. The potential $ is simply 
related to the anisotropy of the CMB on large angular scales via 

5T I 

-7^ - (3) 
T 3 ^ ' 

However, $ is not the most convenient variable for tracing the generation of perturbations 
by quantum fluctuations. For this purpose, it is useful to introduce a second gauge invariant 
quant ity0 

V = a (^ip + (4) 

where 5ip is the perturbation in the scalar infiaton field: ftotaiir, t) = ipQ(t) + 6ip{r, t). During 
inflation, the variable z is 

z = a^/e (5) 
where e is the variable that characterizes the equation-of-state: 

-If^); (6) 



for an infiaton with potential V{ip), pressure p = ^if^ — V{ip), and energy density p = 
+ V{ip). Hence, z = aip/H , where H = d/a = a' /a? is the Hubble constant and 
dot denotes the derivative with respect to the physical time t = J adr. We use everywhere 
the units where AttG = 1. In the post-inflationary phase when the universe is filled with 
hydrodynamical matter, the definition of z in Eq. is replaced hj z = {a/cs)\^ where 
is the speed of sound and, in the definition for v in Eq. (^), 6(p is replaced by the potential 
of the peculiar velocities in the matter. 

By expanding v{x, r) in Fourier modes with comoving wavenumbers k, the equation 
-of- motion for Fourier component Vk becomes (after lengthy computation): 
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{k^ is replaced by c^fc^ in the post-inflationary stage). The canonical variable v used in 
deriving quantum fluctuations beginning from an action describing the scalar fleld coupled 
to Einstein gravity;^ and, it can also be simply related to the gravitational potential $ via 
the constraint equation derived from the — component of the Einstein equations: 



The rigorous method to compute the perturbation spectrum is to solve the second order 
equation for each f ^ using Eq. beginning from when the given wavelength is small com- 
pared to the horizon to when it grows much larger than the horizon. To characterize the 
perturbations we use the power spectrum oi(^= vjz deflned as: 

This power spectrum can be easily related with the power spectrum of the gravitational 
potential if we use Eq. (^), assuming z does not vanish. Note that all of the equations 
above describe the perturbations not only during inflation, but also in the post-inflationary 
(hydrodynamical) Universe. After inflation, e = 0{\) and for the long wavelengths pertur- 
bations (c^/c^ << $fc ~ 0{X)C,k. Hence, is, up to a constant of order unity, the 
power spectrum of the gravitational potential that is of interest for computing perturbations 
of the CMB and the formation of large-scale structure. In this paper, our focus is on the 
fluctuations in and their sensitivity to the equation-of-state, e, during the inflationary 
epoch. 

We note that some discussions of the exact approach improperly characterize C as a 
"conserved quantity." Indeed, an approximate conservation law can be derived for C, using 
Eq. and taking the long-wavelength limit A; — » 0. However the "conservation law" is not 
strictly true for flnite k and neglecting that fact can lead to some confusionJll For example, 
integrating the equation. 



^.c^£VZl±l,,_ (10) 

which follows from — i Einstein equations, one can obtain $ for a given (; in the long- 
wavelength limit {i.e., to lowest order in k), one appears to obtain non-physical, extra 
constants of integration. To see that there are not additional integration constants, one 
should use constraint equation Eq. and keep terms to order k"^. The "extra" constants 
are then fixed in terms of the "old" ones. 

Time-delay Formalism: The time-delay formalism! is one of the methods originally intro- 
duced to compute the energy density perturbation spectrum at the end of inflation. In this 
method, the perturbations are related to the proper time delay 6t{x) between when inflation 
ends at position x compared to the spatial average. Here, we wish to show that this method 
is limited to cases in which the infiaton potential can be treated as nearly flat. To be sure, 
the time-delay formalism is a more intuitive derivation than the gauge-invariant approach 
and treating the potential as flat suffices as a crude estimate for some models. However, if 
one is interested in a rigorous treatment or more general models, including typical models 
of chaotic, natural, and power-law inflation, the time-delay formalism fails. 

The time-delay method was originally presented for a toy model in which the inflationary 
phase is well- described by the de Sitter solution terminated by an instantaneous transition 
to a hydrodynamical stage. The aspects of the toy model which are essential to the time- 
delay approach are: (a) the metric perturbations can be completely ignored during the 
strictly de Sitter stage, and (b) the time delay for the perturbations with comoving wave 
number k, defined as 6Tk = 6(pk{t)/ipo{t), approaches a time-independent constant when the 
perturbation stretches well beyond the horizon. 

Formally, the time-delay formalism cannot be made rigorous for any model, as it follow 
from the equation for the background H = — ^o- The strict de Sitter limit, H = 0, requires 
that ipo = 0, in which case 6Tk = 6ipk(t) / i^o(t) is divergent. However, there is a limited 
range of models for which the formalism gives the leading order contribution. This range is 
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what we wish to clarify here. If one assumes the slow-roll approximation in which (p can be 
ignored in the equation-of-motion for then 

0o~-^,^/3i/, (11) 

where V^^ = dV/ dip. The perturbation 5(f), satisfies the perturbed, linearized Klein-Gordon 
equation: 

6^ + 3H6^ - ^V' + a^o{B - aE)) + V,^^ 6^ + 2V,^ - 0o(0 + 3^) = 0, (12) 

which is supplemented by the — i Einstein equation 

ij + H(t) = ip^Sif. (13) 

In the long-wavelength and slow-roll limits, (p and time- and spatial-derivatives of the metric 
parameters can be dropped: 

msp + {v,^^ -^^y^) = 0. (14) 

(In dropping these terms, we ignore the decaying modes and assume a generic gauge in which 
none of the variables characterizing the perturbations are suppressed compared to others by 
the gauge choice.) The solution for 5ip 

H = c^-f. (15) 

where the constant C = {HV,^ /^)o can be expressed in terms of V and H evaluated at 
horizon-crossing for the given mode; this choice of C guarantees that 6ip H at horizon- 
crossing as expected in the de Sitter limit. 

The key point is that, from Eqs. (|rT]) and (|T^), we find that 

5r = ^ = oc ^. (16) 

Unless V is independent of ipo{t), St depends on time which is inconsistent with an essential 
criterion of the time-delay formalism. Since V is always (/^-dependent in practice, the time- 
delay formalism is, at best, a lowest order approximation. Even so, it should only be applied 
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if the V is extremely flat and H is nearly constant over the range of e-folds of physical 
interest, typically the last 60 e-folds of inflation. If H is nearly constant over the 60-folds, 
then the integral of e = d{l/H)/dt over the last 60 e-folds must be much less than 1/H; if 
the integral is to be less by a factor 5 <^ 1, then we require 

< L (17) 



din e 



dN — 60 ' 

dine ^ 5 . 
^ dN — 3600 ' • • ■ 

where e must satisfy these constraints during the last 60 e-folds, N is the number of e-folds, 
and . . . refers to analogous constraints on higher order derivatives. This represents a narrow 
set of models which excludes common power-law and chaotic inflationary models. 
Horizon-crossing/Bessel Approximation: The horizon-crossing/Bessel approximation 
is based on the exact gauge invariant method but circumvents mode-by-mode integration. 
A recent reviewEl discusses prior work and contains references. The perturbation amplitude 
for a given mode as it enters the horizon in the post-inflationary epoch is expressed in terms 
of the amplitude when the mode crosses beyond the horizon during inflation. To obtain the 
latter amplitude, the solution of Eq. for the non-decaying mode in the long-wavelength 
limit, /c^ < \z"/z\, 

Vk ^ C{k)z (18) 
is matched to the solution in the short wavelength limit, k'^ ^ \z"/z\, 



at horizon-crossing, k = 0{l)aH. (The normalization in Eq. ([T9|) follows from the fact that 
V is the quantum canonical variable and Cg = 1.) The matching condition determines C{k); 
namely, \C{k)\ = 0{1)/V2kz where z = a^p/H evaluated when k = 0{l)aH; hence, the 
power spectrum is 

n(^)-|^|C^(^)P= (l^)^xO(l) (20) 



in the long- wavelength limit, where the subscript means that H and (p are evaluated at 
k = aH precisely. The C(l) ambiguity reflects the fact that the actual matching condition 
is not so precise. 

The Bessel approximation is an improvement of the horizon-crossing approximation in- 
tended to resolve the ambiguity by replacing simple matching at A; ~ aH with a Bessel 
function approximant to Eq. (|^ valid in a range of wavelengths around k'^ = \z" j z\ and then 
matching this solution to the long- and short-wavelength limits. More accurately capturing 
the behavior of the exact solution near k^ = \z" j z\ is important to our purpose because it is 
the integration of Eq. over this wavelength regime that is sensitive to the time-variation 
of the equation-of-state, e(t). Hence, the Bessel approximation not only replaces the 
in Eq. ( [20| ) with a known function, but specifically a function of e. 

The horizon-crossing/Bessel approximation appears at first glance to be a leading order 
expression in a systematic expansion that can be extended to arbitrarily high accuracy. 
Here we show that, in actuality, it is obtained by matching long-, intermediate- and short- 
wavelength solutions. In particular, the horizon-crossing/Bessel approximation assumes that 
the amplitude of a mode is determined by a conditions within a small range of e-folds 
around horizon-crossing for the given mode. If the equation-of-state is changing rapidly, 
this assumption breaks down. 

The ratio z" j z in Eq. can be re-expressed in terms of e: 

z" „^^n 9/ e 3(ilne 1 rflne 1 ( d\v^e\^ l(i^lne\ 

- = 2ffV(^l---j— + + -(^— j (21) 

where < e(A^) < 1 during inflation. Here the conformal time variable r has been replaced 
with A^, the remaining number of e-folds until the end of inflation: the differential dN satisfies 
dN = —aHdr . In Eq. (|^), the cross-over between short- wavelength and long- wavelength 
behavior occurs cit k ^ ~ \z" / z\. Note that this corresponds to horizon-crossing, k ~ aH, 
only provided that e(A^) does not change too rapidly, e.g., {dlne/dN, d'^lne/dN'^) < 0{1)). 
This is a necessary but insufficient condition for the horizon-crossing approximation to be 
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valid. 

Assuming this condition is satisfied, the solution to Eq. for a given mode k can be 
expressed as a Bessel function about k = aH. To see this, it is useful to replace with 

X = \n{Ha/k) = \n{Xph/H-^) (22) 

where Xph = a/k is the physical wavelength of the perturbation with comoving scale k. The 
variable x is roughly the number of e-folds after the given mode crosses the horizon during 
inflation; it is equal to zero at the moment of horizon-crossing, positive after horizon-crossing, 
and negative before horizon-crossing. Then, if we replace v with 

v= {l-ey/^{exp{x/2))v (23) 

Eq. (|^ takes the form 



2 

exp(-2rr) _ 1 (3::::e\'^ _ 3 ding , 1 dln(l-<:) _ 1 ( dlne{l-e) \ _ 1 ln£(l-£) 
4 Vl-e/ 2 dx ^ 2 dx 4 I dx ) 2 dx'^ 



v= (24) 



If e = constant then this equation is reduced to a Bessel equation with exact solutions 
given by V ^ HI ' f j where 

(25) 



2 Vl - 

However, e = constant is not realistic for cosmology since the equation-of-state must change 
near the end of inflation in order to return to the ordinary, Friedmann-Robertson- Walker 
expansion rate. The standard approach has been to solve Eq. (p^ ) approximately near the 
horizon crossing point x = where 

e(x) = eo + eo (^^^) ^ + - (26) 

The subscript is used to symbolize horizon-crossing for the given mode, which is equiv- 
alent to evaluating at x = 0. The solution will be an approximation in which the small 
parameters are eo and x -derivatives, {d.../dx)o. More precisely. In cq is treated as a zeroth- 
order quantity; eo and {dliae/dx)^— are first order quantities; e^, eo {dliae/dx)o, {dliae/dx)^ 
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and {(f\iae/dx'^)o are second order quantities; etc. Note that d\ia{l — e)/dx 
e)){d\ne/dx) is second order. 

Substituting (EBf) in Eq. we obtain: 



d5 

dx 



2 + 



exp(— 22:) 
(l-^o)^ 



1 / 

4 V 1-eo 



(^)o + ^2(x) 



3 /'din 
2 



t;= 



(27) 



where R2{x) denotes the second and higher order terms, 
R2 = 



2e-2- , 3-€o 
T?''' "T 7"; ttJ' 



dln(l-£) 
dx 



1 /dlne\2 1 , f d'^ln 



If dlne \ " in 

4 /O 2 V-"- ' ^-^1 \ dx^ 



3x) 



+ 0{el...) 



(28) 



The initial conditions are fixed by the assumption that the short-wavelength behavior should 
be as in a fiat Minkowski vacuum with only positive frequencies, resulting in the solution 



(29) 



where ^ = exp(— a;)/(l — eo), and 



^~^oY , a f dine 
+ 6 I — — 

dx 



1/2 



(30) 



provided we treat R2 as constant - this is essential to using a Bessel function to represent 
the solution. However, R2 is constant only in the case when e = constant, in which case R2 
itself is precisely zero. The fact that e and, hence, R2 are x -dependent for realistic models 
is what imposes a limit on the validity of Bessel approximation and leads to the result that 
Bessel solution deviates more and more from the exact solution as grows. For example, 
it tells us that it is inappropriate to keep the terms smaller than R2 in the expression for 
index u. Consequently, keeping terms of order eg in the expression for the index u is invalid 
since these terms are comparable to or smaller than the contributions of the x -dependent 
terms in R2 that have been neglected, unless R2 vanishes. Eq. (^0|) should be re-expressed 



3 1 fd\ne\ 



(31) 



where 
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5,vr^\el-\R,{x). (32) 

characterizes the "accuracy" with which the indexes in Bessel function solution should be 
trusted. 

Our analysis illustrates a key limitation to the horizon-crossing and Bessel approxima- 
tions: they cannot be improved further than first order by any reasonable scheme. We have 
just argued that extending the approximation to order requires that we also keep the x- 
dependent terms in i?2; but then Eq. ( |5DD is no longer solvable in terms in Bessel functions. 
Alternative schemes that we can imagine are as difficult as solving the exact equations, and, 
hence, have no advantage. 

The error in ignoring 821' translates into an error in the Bessel approximation to the 
power spectrum that we wish to estimate. To obtain the spectrum, we extend the Bessel 
function approximant valid about x = to the long-wavelength (x ^ or ^ ^ 0) limit: 

^u,(5).£M(l,)-''(i__i_,0(a) (33, 

and the short-wavelength (x <^ or ,^ 00) limit: 

H^;\i) = ^^(p^ + g^)exp - Y - i + (I))) (34) 



where 



^ - ^ 2!{8ip ^ ^yi^) 



(35) 



Combining Eqs. (0), (|31|) and (pSf ) and using H = Hq{1 — eqx + ...), the extrapolation 



of the Bessel solution in the long-wavelength limit becomes 



X ( 1 - ... "^P( + 0(e-3^) ) (1 + 0{52u)) 



4(l-z/o)(l-eo)' 

where j3 = In 2 + 7 — 1 = 0.27 and 7 is Euler's constant, and is a constant to be 
determined by matching to the short-wavelength solution (to be discussed below). The last 
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term in Eq. (|36D characterizes the deviation of the Bessel solution from the exact solution for 
V based on the exact equation, Eq. (^41). This solution to the approximant (Bessel) equation 
is to be matched to the long-wavelength limit of the exact equation {v = C{k)z) to obtain 
C{k) in terms of B^, eg and [dhie/ dx)^. We see that the match is best if x can be chosen 
so that the last two correction factors in parentheses in Eq. (RO) are negligible compared to 



eo and [d\ne/ dx)^. That requires that x be neither too small nor too big where the 
two solutions are matched. The first correction factor requires that the match-point x = Xj^ 
satisfy x+ > max {|ln(eo; (^)o|}' second factor requires that a;+ be sufficiently 

small that 

52y{x+) « min(eo; {^-^^ )• (37) 

(Recall that 82^ includes x-dependent terms that increase in magnitude with increasing x.) 
The result is that 



~ max 



In 



mm(eo; 



dine 

dx 



OJ 



(38) 



is optimal for obtaining the best match. The fact that the optimal match-point, x+, is 
constrained from above and below means that there is a residual second-order error which 
cannot be improved upon by this matching procedure. 

The matching of the Bessel solution to the short-wavelength limit proceeds similarly and 
determines Bk- Eq. (|^) implies that 



V = Bus^ - + 0{l) ^^^J a; + -exp(2a;) + 0(e"")Jexp(z(fcr). (39) 

where the second terms inside the brackets comes as a result of uncertainty in Bessel function 
index, 821' ■ This solution can be matched to the short-wavelength limit, v ^/T/2k exp{ikT) 



, roughly 



to first-order accuracy in eo and {d\iae/dx)o at match-point x„ ~ ^In ^ '^^'^^j^T^"^ 
two or so e-folds before the end of inflation for typical models. (Note that Eq. (^) con- 
tains two second-order correction terms, one of which increases with x and the other which 
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decreases with x. Consequently, just as with we find that the point x_, the short- 
wavelength match-point which gives the best accuracy, is constrained from above and below.) 
The higher-order corrections are negligible provided 



6oB(x) = 0(1] 



dln{l — eo) 
dx 



In 



(iln(l — eo) 
dx 



< 1. 



(40) 



In this limit, = A/vr/4/c(l + 0{l)62B(x^). Substituting this expression for Bf^ into 
the long- wavelength expression Eq. ( p6D and obtaining a complete first-order expression for 
C{k), one can then determine the power spectrum: 

dine' 







2 






27r2 


z 









2/?e+(l-/3) 



dx 



0{l)62B{x_) + 0{l)62iyix+ 

(41) 



where the last two terms in square brackets characterize the uncertainty and should be 
smaller than the previous, first-order terms. Comparing to the original estimate, Eq. (pO|), 
we find that the 0{1) factor has been replaced by a function which depends on the equation- 
of-state. The spectral index of the scalar power spectrum is then 

dlnPr 



ns = l + 



dink 



dine n ^ dine , ^.(P'lne 
1 - 2e - -— - - 2/5e-— + 1 - /? 

dx dx dx^ 



+ . . . 



(42) 



and 

dris 
dx 



-2e- 



d{ns - 1) 
dx 

dine (ff'lne 



dx dx^ 



Ae' 



,dlne 
dx 



2(3e 



dine 

dx 



2(3e 



rf2 Ine 

dx"^ 



+ (1-/5) 



d^ ln( 

dx'^ 



(43) 
+ ... (44) 



where . . . are uncertainty due to and 52B. These expressions agree with previous 
result sIll'llZl except that they are expressed in terms of the equation-of-state and its deriva- 
tives. 

Eqs. ( P?!) through ( PI ) summarizes the basic result: the Bessel approximation for P(^{k) 
is only good to first order in eo and {dlneo/ dx)^ at best, and then only if 52B{x^) and 
52Z/(x+) are negligible compared to the first order contributions. In previous discussions, it 
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was pointed out that the approximation was only vahd if e and d\n.e/dx (or equivalents) 
are nearly constant over some range of e-folds around horizon-crossing but, otherwise, 
the conditions for the approximation to be valid were not specified. Here, we see that the 
relevant range of e-folds is between x_ and x+ e-folds about k = aH, typically the five or so 
e-foldings surrounding horizon-crossing, k = aH. We also see that, clS X ^ X-^ , the Bessel 
solution approaches the exact solution to within accuracy 621^ and 62B] for x beyond this 
range, the Bessel approximant diverges from the true short- and long-wavelength solutions. 
Consequently, the Bessel approximant can achieve first order accuracy, but no better. Recall 
that we must also restrict ourselves to d"^ In e / dx"^ < 0{1) (see discussion under Eq. (^I]) . 

In special satisfactory numerical result can be obtained even though some con- 

straints are not satisfied: A prominent example is natural inflation and other potentials of 
the form V ^ Vq — aip'^ + . . .. In Eq. (|38D, two independent constraints are implied by the 
parenthetical (eo; (^^)o)- '^^^ ^^^^ constraint ensures that the Bessel approximation gives 
the correct result to leading order in eg; the second ensures that the correct result to leading 
order in (^j^. For some natural inflation models, though, the first constraint is strongly 
violated so that the O(eo) terms in the Bessel approximation cannot be "trusted." However, 
not only is the second constraint satisfied, but the terms in the Bessel approximation propor- 
tional to {dlne/dx)o are so much larger than the O{eo) terms that the violation of the first 
constraint is numerically insignificant. The success of the Bessel approximation is accidental 
in this sense (and may have deceived some into thinking that the Bessel approximation has 
a much wider domain of validity than it actually does). 

Our elaborate analysis can be reduced to a simple statement about domain of validity 
of the horizon-crossing/Bessel approximation: 

(1) If e = constant, any value in the inflationary range between and 1, the Bessel solution 
is exact. However, a model with e = constant is not physically realistic since inflation never 
terminates. 
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(2) If e 7^ constant, the Bessel approximation is only accurate if eo and {d\ne/dx)o are small 
enough that the first order contributions in eo and {d\ne/dx)o are much larger than the 
higher order contributions, 621^ and 62B. Suppose we demand that the higher order terms 
on left-hand-side of Eq. (^) be less than a factor 5^1 times the first order terms on 
the right-hand-side of Eq. (^Sp. (For example, 6 might be determined by the resolution of 
an experiment and we may wish to know if the Bessel approximation provides the needed 
accuracy.) Assuming no accidental cancellations, Eqs. (|37|) through (|40|) reduce to: 



where x+ = max ||ln(eo; (^)ol} > ^ ■■■ refers to analogous constraints on higher 
order derivatives. Even for modest accuracy, 6 = 20% and x+ ~ 2, the ratio is ~ 0.1, 
enough to highly restrict the range of e and its derivatives. 

Comparing the constraints above to Eq. (|l^, one sees that the horizon-crossing/Bessel 
approximation applies to a wider range of models than the time-delay formalism. Never- 
theless, the range is narrow compared to full spectrum of inflationary models. One class of 
models in which the horizon-crossing/Bessel approximation is valid, where {dlne/dx)o -C 
eo < 5/x+ <^ 1, includes the simplest models of new infiation&iJii, chaotic infiationlii with (f)"' 
potentials and n >> 2, and extended inflation0i3, which are realistic models incorporating 
inflation. For these models, one obtains the CMB anisotropy prediction: and r obey the 
relationii: r ~ 21(1 + 7) ~ 7(1 — Us), where r ^ e/14 is the ratio of the tensor mode to 
the scalar mode in terms of the contribution to the CMB dipole moment. The second class, 
where eo -C {dlne/dx)o ^ ^ 1, includes a range of natural inflation modelsj^l chaoti 
inflation models with 0" potentials and small n, and some two-fleld inflation modelsB i 
which the inflaton fleld rolls near an extremum of the potential during inflation. 

Some may have assumed that the good agreement between the Bessel approximation and 



ic 



m 
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the exact methods for these two cases meant that the Bessel approximation could be used 
for a broader range of models. In fact, our results show that these are essentially the only 
models for which the approximation can be trusted. 

Impact on Microwave Background Anisotropy Prediction: An Illustration The 

error in using one of the approximate procedures instead of the more cumbersome mode-by- 
mode integration propagates to predictions of the cosmic microwave background anisotropy 
and large-scale structure. As a dramatic illustration, Figure 1 shows a comparison of the 
predicted CMB anisotropy power spectrum using the time-delay formalism or naive horizon- 
crossing approximation (based on Eq. {^^), the Bessel approximation (based on Eq. (^l]) ), 
and the exact computation for a sample inflaton potential, V{(f)) = A^(l — | tan^^ {5(f> / rrip) , in 
which the equation-of-state changes rapidly enough near = that Eq. (^Sj) is not satisfied. 
(For this toy model, we have taken (p ^ —0.3 to correspond to 60 e-folds before the end of 
inflation.) The discrepancy in the CMB predictions is large compared to the anticipated 
experimental resolution of future space-based anisotropy experiments. Less dramatic effects 
occur in more typical models with slowly varying equation-of-state; an analysis for a wide 
spectrum of models be presented in a future paper. @ 

Summary: Our conclusions are summarized in Eqs. (0) and (^5|) as constraints on the 
equation-of-state, e. The basic result is that the time-delay and horizon-crossing methods are 
reliable approximations only if e and its time-variation are rather small. These constraints 
can be re-formulated in terms of rules-of-thumb for an inflaton potential, V{(f): Assuming 
higher-order corrections to our approximation should he 6 < 0.20% and x+ ~ 2, then, if 
V{ip) satisfies any of the following conditions (recall that AttG = 1 and > 2): 




(47) 



(46) 



(48) 
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during the last 60 e-folds of inflation, the horizon-crossing/Bessel approximation is not 
rehable and mode-by-mode integration is required. For the time-delay formalism, the con- 
straints are roughly 60 times more stringent. 

An important consequence is that attempts at precise fltting of CMB anisotropy 
data@~0 and large-scale structure measurements and attempts to "reconstruct" the in- 
flaton potential from CMB, as described in a recent review,^ is not as straightforward 
as one hoped. If the horizon-crossing and reconstruction approaches were generally valid, 
then inflationary predictions could be parameterized with only a few variables ( e.g., cq and 
{dlne/dx)Q evaluated for the mode crossing the horizon in the present epoch). Simultaneous 
fitting of these parameters along with other cosmic parameters, (such as the Hubble con- 
stant, the cosmological constant, the baryon density, etc.) would provide tight constraints 
on all. Indeed, this approach has been assumed in most prior discussions of fltting data. 
To be sure, cases where the equation-of-state is nearly constant and the horizon-crossing 
approximation is valid appear to be the simplest forms of inflaton potential based on our 
current understanding. So, one can decide a priori to assume this subclass of inflationary 
potentials; in this case, there is no point to general reconstruction methods since the po- 
tential forms are set by the a priori assumption. Alternatively, one may make no a priori 
assumptions, in which case reconstruction methods are not useful since they are not valid for 
general potentials. If we broaden the spectrum of possible potentials, the fltting of cosmic 
parameters must be learned by comparing data to some systematic search through exact 
results obtained by mode-by-mode integration. How best to perform the search and how this 
affects the empirical resolution of cosmic parameters from CMB measurements is a subject 
of current investigation.!! 

This research was supported by the Department of Energy at Penn, DE-FG02-95ER40893 
(LW and PJS), and by the Tomalla Foundation (VM). 
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FIG. 1. A comparison of the horizon-crossing and Bessel approximations to exact 
mode-by-mode integration for an inflaton potential in which the equation-of-state (e) is varying 
rapidly. The power spectrum has been computed and converted into a prediction of the CMB 
temperature anisotropy spectrum on large angular scales. 
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